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ABSTRACT 

We analyze the production in e"^e~ collisions of new heavy fermions stemming 
from extensions of the Standard Model. We write down the most general 
expression for the production of two heavy fermions and their subsequent de- 
cays, allowing for the polarization of the e"'"e~ initial state and taking into 
account the final polarization of the fermions. We then discuss the various 
decay modes including cascade and three body decays, and the production 
mechanisms, both pair production and single production in association with 
ordinary fermions. 
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1. Introduction 



Despite its tremendous sucess in describing all experimental data available today, the Stan- 
dard Model of the electromagnetic, weak and strong interactions based on the gauge sym- 
metry SU(3)c X SU(2)l X U(1)y, is widely believed not to be the ultimate truth. Besides 
the fact that it has too many parameters which are merely incorporated by hand, the 
Standard Model does not unify the electroweak and strong forces in a satisfactory way since 
the coupling constants of these interactions are different and seem to be independent. There- 
fore, one expects the existence of a more fundamental theory which describes the three forces 
within the context of a single gauge group and hence, with only one coupling constant. This 
grand unified theory will be based on a gauge group containing SU(3)c x SU(2)l x U(1)y 
as a subgroup and will reduce to this symmetry at present energies. 

The grand unified groups [1-4] provide fermion representations in which a complete gener- 
ation of Standard Model quarks and leptons can be naturally embedded. In most of the 
cases these representations are large enough to accomodate additional new fermions which, 
in fact, are needed to have anomaly-free theories. It is conceivable that these new fermions, 
if for instance they are protected by some symmetry, acquire masses not much larger than 
the Fermi scale. This is very likely, and even necessary if the new gauge bosons which are 
generic predictions of the unified theories, are relatively light |p. 

Besides the SU(5) group [the simplest Lie group containing SU(3)c x SU(2)l x U(1)y as 
a subgroup and with two representations to accomodate the 15 Standard Model fermions] 
which has no room for "light" new fermions or gauge bosons, the SO(IO) group ^ has 
received much attention. It is the simplest group in which the 15 Weyl spinors of each 
Standard Model generation of fermions can be embedded into a single multiplet. This rep- 
resentation has dimension 16 and, to have an anomaly-free theory, contains a right-handed 
Majorana neutrino. In fact, heavy isosinglet neutrinos have been discussed in various mod- 
els, such as left-right symmetric models [^, in attempts to explain the small masses of the 
three observed neutrinos. 

Another popular unifying group is Eg [^ which contains SU(5) and SO(IO) as subgroups 
and is the next anomaly-free choice after SO (10). The interest in Eg is mainly due the fact 
that superstring theories, which attempt to unify all fundamental forces including gravity, 
suggest that this symmetry is an acceptable four dimensional field theoretical limit In 
Eg, each quark-lepton generation lies in the representation of dimension 27; to complete this 
representation, twelve new fields are needed in addition the Standard Model fermion fields. 
For each family one has two additional isodoublets of leptons, two isosinglets neutrinos [which 
can be either of the Dirac or Majorana type] and an isosinglet quark with charge —1/3. 

Several other gauge groups have been considered with various theoretical motivations and 
most of them predict the existence of new fermions. For instance, schemes of grand unifi- 
cation based on large orthogonal groups have been proposed to explain the origin of parity 
violation in weak interactions [^: in these models, weak interactions are parity symmetric 
but fermions with left-handed and right-handed couplings aquire different masses. They 
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predict a rich spectrum of fermions, called mirror fermions ||^ , which have the opposite chi- 
ral properties of the ordinary ones. In the simplest version of these models [1^, the gauge 



symmetry and the symmetry breaking pattern are the same as in the Standard Model: one 
simply adds to the spectrum of the latter three families of heavy fermions with opposite chi- 



ralities. Theoretical arguments based on the unitarity of scattering amplitudes suggest 
that the masses of these mirror fermions should not exceed a few hundred GeV. 

The direct search for these new fermions and, in case of discovery, the study of their basic 
properties, will be a major goal of the next generation of accelerators. In this paper, we 
analyze in detail the production of these new fermions at e^e~ colliders. 

The new leptons and quarks will mix with the ordinary fermions of the Standard Model 
|1T2| , |T^ . This mixing will give rise to new currents which determine to a large extent their 
decay properties and allows for new production mechanisms. If the new particles have non- 
zero electromagnetic and weak charges, they can be pair produced if their masses are smaller 
than the beam energy. In general the reactions are built-up by a superposition of photon 
and Z boson exchange [additional contributions could come from extra gauge bosons if their 
masses are not much larger than the cm. energy of the collider]. The cross sections are large 
H and, up to phase space suppression factors, of the order of the point-like QED cross 



section for muon pair production. 

Fermion mixing allows an additional production mechanism for the new fermions: single 
production in association with their light partners. In the case of quarks [and for second and 
third generation new leptons if inter-generational mixing is neglected] the production process 
is mediated by s-channel gauge boson exchange; since the mixing angles are restricted to 



values smaller than 0{1Q ) by present experimental data |T^, the cross sections are rather 
small. But in the case of [the first generation] heavy leptons, additional t-channel exchanges, 
W exchange for neutral leptons and Z exchange for charged leptons, are present increasing 



the cross sections by several orders of magnitude. This results in large production rates [ p^ | 
which permit to probe lepton masses close to the total cm. energy. 

The new fermions will decay through mixing into light fermions and gauge bosons. De- 
pending on the particle masses, the gauge bosons can be real or virtual and will decay into 
light quarks and leptons. Therefore, in the production of the new fermions the final states 
are rather complicated: six particles in the case of pair production and four particles for 
single production. However, it is very important to have at hand the information on the 
correlation between all the particles involved in the process. Indeed, these correlations will 
be very helpful to optimize the experimental cuts which permit to suppress the various back- 
grounds without affecting drastically the signal cross sections. Furthermore, they permit to 
discriminate between particles of different nature [e.g. Majorana or Dirac neutrinos] or with 
different couplings [vector, "mirror" or standard couplings in pair production, or with left 
and right-handed mixing in single production] and therefore shed some light on the origin 
of the new fermion. 
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Several analyses of the production of new fermions in e^e" collisions have been conducted 
in the recent years [14-20] in various special cases. In this paper we will extend on these 
analyses in the following ways: 

{i) We present the most general expressions for the production of two heavy fermions, with 
different masses to allow for the production of two different fermions and to treat the cases 
of single and pair production in the same footing, of any flavor and for arbitrary couplings 
of the new fermions, including all possible channels and the polarization of the initial e^e" 
state; complete formulae are given for angular distributions and total cross sections. 

{ii) We give the most general expression of the decay of the heavy fermions into three body 
final states, allowing for cascade decays, including all channels and the possibility of off-shell 
intermediate vector bosons; complete formulae for the angular and energy distributions of 
the final decay products as well as the total decay widths of the heavy fermions are given. 

{Hi) We systematically take into account the polarization of the heavy fermions in both 
the production and the decay processes; due to the factorization of the two sequences, the 
information on the final polarization permits an easy reconstruction of the full correlations 
between the initial state and the final particles from the decays of the heavy fermions. 

Note that since our analytical results are general, they can also be used to discuss, at e^e~ 
colliders, the details of top quark pair production, or the production of heavy fermions of a 
fourth generation with a heavy neutrino. 

The paper is organized as follows. In the next section, we summarize the interactions of 
the new fermions. In section 3, we analyze the general case of the production of two heavy 
fermions and discuss the cases of pair production and the single production in association 
with ordinary fermions. In section 4, we discuss in detail the decay modes of the new fermions 
including the cascade decays and the three body decays with off-shell gauge bosons. Section 
5 contains our conclusions. For completeness, we summarize in the Appendix the formalism 
for combining the spin-dependent cross sections and decay distributions which permits to 
obtain the full correlations between all the particles involved in the process. 

2. Interactions 

The new fermions couple to the photon and the clectroweak gauge bosons W/Z with full 
strength, except for singlet neutrinos which have zero electromagnetic and weak charges and 
therefore couple to the latter only through mixing as will be discussed later. They similarly 
couple to extra gauge bosons when the latter are present. Allowing for these extra currents, 
the interaction is given by the Lagrangian 

^= E 9vJ^V^ (2.1) 
v=^,z,w,- 
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where the complex conjugation of the charged currents is understood. The currents can 
be expressed in terms of left-handed and right-handed charges as 



f 



Qff'V (1 - 75) ^ Qff'V (1+75) 



(2.2) 



2 "'^ 2 

For the minimal Standard Model gauge bosons, the coupling constants gy are simply 



g^ = e = V Area , gz = e/swcw , gw = ev2/sw (2.3) 
with e the proton charge and = 1 — = sin^ 6w', the coupings Q{Jji read 

^L,R - e , (^L,R - ^3L,3R " ^ , ^l,R " MsL.SRI 

with e-^ the electric charge in units of e and I3L, I3R the third components of weak isospin. 

The inclusion of additional gauge bosons in the previous equations is straightforward once 
their couplings to fermions are specified. 

The new fermions will mix with the ordinary fermions which have the same U(1)q and SU(3)c 
quantum numbers [1^, [13[. The mixing will determine the decay systematics of the new 
leptons and quarks and allows for their single production in association with light fermions. 
In principle, one has to treat the three generations of fermions on the same footing; this, 
leads to rather complicated mixing patterns. For instance in Eg, the mixing in the general 
case would be described by 6 x 6 non-diagonal matrices in the quark and charged lepton 
sectors and by 15 x 15 matrices in the neutrino sector |^. The matrix elements will depend 
on the vacuum expectation values of the Higgs fields and on arbitrary Yukawa couplings, 
making a general analysis rather complicated^. 

In order not to commit ourselves to any particular model, we will allow for the mixing 
between different generations and treat the mixing angles as phenomenological parameters. 
For instance, in the case of the interaction of the electron and its associated neutrino with 
exotic charged and neutral heavy leptons, the general Lagrangian describing the transitions 
between e, z/g and the heavy leptons A^'^, where is a generation index reads 
3 

2 



k=l i=L,R 
3 



+ h.c. 



+ \T. [aw Cr^'^e7,iVf + gz Cfi^el.N^ Z^] + h.c. (2.5) 



2 

^ fc=l i=L,R 



where we have allowed for both left-handed and right-handed mixing and assumed small 
angles so that one can write sin ~ Cl,r- From this Lagrangian, where gw,z are given in 
eq. (2.3), the charges Q{^/ as in eq. (2.2) can be easily derived. The generalization to the 
other light leptons and to quarks as well as to extra gauge bosons is straightforward. 



^However, the inter-generational mixing will induce at the tree level, flavor changing neutral currents 
which are severely constrained by existing data 0| ; neglecting the latter allows for an enormous simplification: 
the mixing can be parametrized by a few angles which can be treated as phenomenological parameters. 
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3. Production in e+e Collisions 

3.1 General Case 

In this subsection, we give the most general expression of the differential cross section for 
the production of two fermions with different masses, to treat pair and single production on 
the same footing, including the longitudinal polarization of the initial je^ beams and the 
polarization of the final fermions. 

Consider the process where a pair of heavy fermions is produced in e^e~ annihilation through 
gauge boson exchange 

e+(7,e) e-(/,0 ^ F(p,n) (3.1) 

^, ^ denote the degrees of longitudinal polarization of the initial electron and positron and 
/, I their momenta; n{n) and m{fn) are the four-momentum, spin vector and mass of 

the final fermion F{F). The most general form of the differential cross section, da can be 
written in terms of the scalar products of the spin four-vectors and momenta of the particles, 
as 

where e is the proton charge, Nc the color factor of the final fermions and the squared 
amplitudes A and A' can be expressed in terms of generalized charges Qi,2,3 and Q[ 2 3 

A = {pri){pd){Q, + Q,) + {pd){pJ){Q,-Q,)+mm{iri)Q2 

- in.l)[mip.l)Q', + m(p.I)(g; - Q'3)] + in.l)[m{p.l)Q', + m{p.l){Q\ + Q'3)] 

- {n.l)[m{pri)Q'^ + m{p:i){Q\ + Q'^)] + {nri)[m{p.l)Q'^ + m{p.l){Q\ - Q'^)] 

+ {nn)Q2 [(/.I) {p.p) - {p.l) (p.I) - {pi) (p./)] - 2Q2 [{n.l) (n.l) + {n.l) {n.l)]{l.l) 

- {n.l)(n.l)[Q2{{Ll) + {p.p) + {p.l) + (p.I)) + mM{Q^ + Q^)] 

- {n.l){n.l)[Q2{{l.l) + (p.p) + (p.l) + {p.l)) + mm{Qi - Q,)] 



A' = A{Qi^Q[ , Q2^ Q2 , Qs ^ Q3 ) 
In terms of the helicity amplitudes Qij with i,j = L, R, the charges Q and Q' are 

Qi = \[\Qll\'+\Qrr\' + \Qrl\'+\Qlr\ 
Q2 = -R-e [ QllQ*rl + QrrQlr] 
Q3 = 1\\Qll\'+\Qrr\'-\Qrl\'-\Qlr\ 



(3.3) 
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Q2 



\Qll\'^ + IQrl^ — \Qrr\'^ — \Qlr\ 



Re [ QllQ*rl - QrrQlr] 



Q3 = 7 IQllI'+IQ 



LR\ 



\Q 



RR\ 



\Q 



RL\ 



(3.4) 



The hehcity amphtudes Qij depend on the process under consideration. For s — (/ + 1)^, 
t — {p — Vf and u — {p — ly channel vector bosons V5, Vt and Vu exchange, respectively, 
the general form is 



I 9Vt ^Qj ^ ^ I 3Vjj Qi "Qj " ^ 

^ff'v 



(3.5) 



The normalization factors gy and the reduced couplings Qifji can be derived from the 
Lagrangian describing the ff'V interaction; see section 2. 

Integrating over the variables of one of the final fermions as well as on the azimuthal angle 
of the remaining one, the differential cross section dcr/dcos6', where 9 specifies the direction 
of the latter particle with respect to the incoming electron, reads 



dcost^ 8 4 L 



(3.6) 



where ao = iira'^/Ss is the point-like QED cross section for muon pair production and A the 
usual two body phase space function 



s 

+ 2— n- / + / 

s 

tTT — 

+ 2jn-{l-l) 

Wi — 
- 2— n •(/ + /) 
s 



X = {1 — /j,^ — jj^Y ~ ^l^^l^ 1 with jj, — mj -sfs , = mj \fs 
In terms of the charges Qi and the reduced amplitudes squared A and Al read 

A ^ 1 - (^2 _ -2^2 _^ ^Qg2 ^1 _^ 4/i7Z(52 + 2A5 cos BQ'i 

2-n-{l-l) \{1 + ]l')Q[ + (1 + - 7l')-Q^ -\ 
L m 

1 777 

(1 - + ji^)Q'^ + x-2 cos ^(g; - -g^) 
(1 - 71^ + //^)g; + (1 + /i^ - //^) =g^ + cos eg'3 

(1 - /Z^ + ;,2)g'^ + COS ^(g; - -Q'^) 

Tn 

8 

+ n-n (1 -cos^6') A g2 Q2[n-l n ■ {I + I) +n ■ I n ■ {I + I)] 



(3.7) 



+ -n-ln-l [{1 + 1/ + Ji'^ + cos e)Q2 + 2iJ4i{Q3 - Qi) 
s 

4 - 1 

- -n • Z n • / [(3 - /i^ - /Z^ + A5 cos^)g2 + 2^Jl{Q:i + gO 

5 
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A = A{Qi^Q[,Q2^ Q'2 , Qs ^ Q3 ) (3.8) 

The polarization four-vector of the final state fermion F is defined by d(jP°Ydcos0 ~ 
dcr^^'^P°Yd cos 6' X [1 + Pfj^n^], with the spin vector which satisfies the relations n ■ n = —\ 
and n-p = 0; see Appendix. In the F rest frame, assuming CP-conservation, the components 
are (0, Pj_, 0, P||) with Pj_ and P\\ the transverse and longitudinal polarizations with respect 
to the flight direction. Summing over the polarizations of one of the final fermions, e.g. F, 
the longitudinal and transverse components of the polarization vector of the other fermion, 
in its own rest frame, are given by 

[i-/.^ + 7z^ + (i + - Tf) cos2 e]\iQ', + COS + A)|g'2 + 2(1 - /^^ - 7^')Q'l] 

n I = 



1 - (^2 _ ^2)2 + X COs2 0] Qi - 4^/1^2 - 2A^ COS 
■2\mj 



Ti')Q'i + (1 + - 7i')^Q'2 + A2 COS 0Q^j 



[1 - (/x2 - 772)2 + ^ gog2 ^] + 4^7IQ2 + 2A 2 COS OQ^ 

In the next subsections, the special cases of heavy fermion pair production and single pro- 
duction in association with massless fermions will be discussed. 

3.2. Pair Production 

In e+e^ collisions, the pair production of new fermions proceeds through s-channel gauge bo- 
son exchange; there are also contributions from t-channel exchange in the case of heavy lep- 
ton production, but they are quadratically suppressed by mixing angle factors and therefore, 
rather small. The unpolarized differential cross section da/d cos 9 for the process e+e" — >■ FF 
is (see also Ref. [pO(] ) 



1 Q 

- (To AT,/?^ f(l + /5|cos2^)Qi + (l-/3|)Q2 + 2/5i.cos0Q3l (3.10) 



dcos^ 8 

with (3f = — Arn^p/ sY^"^ the velocity of the fermion in the final state; the charges Qi, Q2 
and Qs are given by in eq. (3.4) with the helicity amplitudes Qij with 2,j = L, i? in the 
general case in eq. (3.5). If only s-channel photon and Z boson exchange is present, these 
helicity amplitudes are simply given by 



QFFZQe ^ 



where Q{^^ are the reduced couplings of the left and right-handed fermions to the Z boson. 

In the case of unpolarized initial beams, the cross section eq. (3.10) allows for three indepen- 
dent measurements: the total production cross section ctf-, the forward-backward asymmetry 
A^^ and the parameter ap defined as dcr/dcos^ ~ 1 + a^^cos^^. These three parameters, 
and their asymptotic values, are given by 
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CTF= icjoN^Pp ^i + l[3l)Q, + [l-[3l)Q 



Zip 



(3.12) 
(3.13) 
(3.14) 



For Majorana neutrinos one has to symmetrize eq. (3.10) because of the two identical par- 
ticles in the final state. This symmetrization makes that the Majorana neutrino has only 
axial-vector couplings so that the charges are Q2 = —Qi and = 0, leading to the simple 
expression for the cross section 



\Q 



LL\ 



\Q 



RL\ 



(3.15) 



The total cross section is proportional to and thus, strongly suppressed near threshold; the 
angular distribution behaves like da/dcosO ~ 1 + cos^^ and therefore, there is no forward- 
backward asymmetry and the a parameter is equal to one. Note that since the isosinglet 
Majorana neutrinos do not couple to the photon and Z boson, they can be pair produced 
only through the exchange of an extra gauge boson and thus, the helicity amplitudes Qij 
have to be altered. 

Let us now discuss the polarization of the heavy fermions. Summing over the polarizations 
of -F, the two components of the polarization vector of F in its own rest frame, P\\ and Pi 
can be written as (see also refs. [pOl, p2|) 



1 + cos^ e)(3FQ'^ + cosg [(1 - (3l)Q'^ + (1 + (3l)Q\] 
(1 + 131 cos2 e)Q^ + (1 - Pl)Q2 + 2pF cos ^Qa 



1-I3lsme- 



1 + pl cos2 e)Qi + (1 - [3l)Q2 + 2I3f cos OQ. 



(3.16) 



where the helicity amplitudes Qij are given in eq. (3.5) and (3.10) and the charges <5'i2 3 
eq. (3.4). Averaged over the polar angle, the two components become 



m 



< Pii > 



3(l + /?2/3)g^ + 
37r mi? Q'l 



4 yi(l + /?|/3)gi + (l-/?2)g2 



(3.17) 



Again, for Majorana neutrinos one has to symmetrize the previous expressions. This makes 
that the currents become purely axial- vector, and therefore there is no polarization effect. 
The polarization of the final particles, together with their angular distributions, permits 
to discriminate between Majorana and Dirac neutrinos, or between particles with different 
couplings [standard, mirror or vector couphngs]. 
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4. 2. Single Production 



In e^e~ collisions one can also have access to the new fcrmions via single production in 
association with their light partners if fermion mixing is not too small. Assuming that extra 
gauge bosons are too heavy to affect the production, the process proceeds through s-channel 
Z exchange for all fermions, but for heavy leptons [only the first family if inter-generational 
mixing is neglected] one has additional ^-channel gauge boson exchanges: W exchange for 
neutral leptons and Z exchange for charged leptons. Neglecting the mass of the light fermion 
partner, the differential cross section for the process e+e" Ff is 



da 



dcos^ 8 



CTo iVe(l 



(1 -/x^)cos^^ Qi + 2cos^g3 



(3.18) 



where /i^ = rrip/s. The charges Qi^s as given in eq. (3.5) are built-up by the helicity 
amplitudes 



Qij — 


1 cr^'^of ^ 








1 - 


Qij — 




1 




- z 


Qij — 


1 crQf 


1 




— z 



+ 



1 



+ 



t/s — z 

cN'^WqciAV 



t/s 



w 



for E' 



for N'' 



for quarks 



(3.19) 



with the reduced masses z — M'^/s, w — M^/s and t/s — —(1 — cos^)(l — /^^)/2. Note 
again that in the case of Major ana neutrinos, the Qij have to be symmetrized. 

For quarks [and for second and third family leptons if inter-generational mixing is neglected] 
the total cross section (7(e+e~ — > Ff) takes the simple form 

1 



cr e ' e 



2\2 



Ff) = ao N,{1 - ii') 



1 + 2^ 



(3.20) 



The cross sections for the production of the conjugate states a{e^e'' — > Ff) is the same. For 
first generation heavy leptons, the analytical expressions are much more involved because of 
the t-channel contributions. Denoting by Vs the gauge boson exchanged in the s-channel 
[Z for both types of leptons] and by Vt the one exchanged in the t-channel \W for N and Z 
for E], one can write a common expression for heavy lepton L = N^, E^ single production 



a (e^ 



U) = 3ao -(1 



2\2 



y+ "T <i- 



+ 



+ 



(1-//')(3 + 2^;t-//') -2(1-//' + ^;t)(1 + ^;t) log 



I- [J? +Vt 



Vt 



ysVr 



Vs 



-(1 - At')(l + 2vt) - 2vt{ii'' - vt) log 



1- fl^ +Vt 



Vt 



.VsVt 



Vs 
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+ 



+ 



1 + 2vt 

Vt 



(2 
1 



/i^ + 2vt) log 



1 - /i^ + f T 



Vt 



-VtVt 



1 + f r - /i^ 



- {2vt - /i^) log 



Vt 



^VtVt \ 

<1. j 



(3.21) 



where Vs = My^/s, Vt = My^/s and the charges q^' = gX^^^ ' ' ^^e defined by 



W 
5+ 



I VV'i2 I I VV i2 



VV 



1 r 
4 



I VV'i2 I I VV'i2 



with q^""^ = {glJe')Q^'''^Qf' etc. 

To obtain the cross section for the production of E and in the case of left-handed or 
right-handed mixing is then straightforward: one has simply to specify Vs and Vt and 
choose the proper combination of charges. For instance, in the case of the production of a 
heavy neutrino with a left-handed mixing denotedQiVL, one recovers up to a factor of two, 



the formula given in Ref. |18| for the production of Majorana neutrinos, by setting 



Ml 



Vt 



WW 



wz 







(3.22) 



This is due to the fact that the production cross section of a Majorana neutrino is just the 
sum of the production cross sections of a Dirac neutrino and its anti-neutrino; this also holds 
true for the angular distributions. 

For E production, we have a Z exchange in both s and t-channels. This simplifies consider- 
ably the expression eq. (3.21) since in addition to the fact that vs = vt, the charges factorize. 
Note that for s^^ = 1/4, the left-handed and right-handed couplings of the electron to the 
Z are equal in magnitude but with opposite signs. This leads to 



1- 







(3.23) 



which translates into the fact the cross sections for El and E^i are approximately equal. 

Finally, the longitudinal and transverse components of the polarization vector of the heavy 
fermion produced in association with a massless partner reads 



P|l 
Pi 



[1 - /i^ + (1 + /i2) cos^ 9]Q'^ + 2 cos eQ[ 
[1 + + (1 _ ^2) cos2 e]Qi + 2 cos eQ3 
Q\ + cos OQ'^ 



2/i sin^ 



[1 + /i2 + (1 - ;i2) cos2 e]Qi + 2 cos eQ3 
where the charges Q'123 given in eq. (3.4) with the helicity amplitudes eq. (3.19). 



(3.24) 



^In principle, the indices L, R refer to the handedness of the heavy lepton mixing with its hght partner. 
However, due to the fact that the latter is massless, they are also the chirality of the heavy lepton. 
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4. Decays and Correlations 



4.1 Two-body Decays 



The heavy heavy fermion F, with a mass m and spin four-vector n, will decay into a lighter 
fermion /o and virtual or real gauge bosons which subsequently decay into two massless 
fermions /i and /2 



F{p,n) 



/o(/o) V{lv) 



/o(/o) /i(/i) /2(/2) 



(4.1) 



In most of the cases the fermion /o is just the ordinary partner of F [which can be considered 
as massless] and the decay occurs through mixing. But it is possible that /o is also a heavy 
fermion and the process eq. (4.1) is a "cascade" decay. For instance, for fermions belonging 
to the same isodoublets, the heavier fermion can decay through the exchange of a boson 
into its lighter isospin partner and the latter subsequently decays through mixing into three 
massless particles. Moreover, it is also possible that the mixing between heavy fermions 
is much larger than the mixing between the new and ordinary fermions in which case the 
decay of the heavy fermion first into a lighter one, is more important although kinematically 
disfavored. 

If the mass difference between F and /o is larger than the mass of the exchanged gauge 
boson, the latter will be on mass-shell and the decay is a two-body decay. Assuming that /o 
is also polarized [the mass and spin four-vector will be denoted by mo and no respectively], 
the differential decay width in terms of the momenta and spin-vectors of F and /o is given 

by 



dr= (27r)^5^(p-Zo-/ 



dHo dHv 



V 



(27r)32/o (27r)32/o 2m 



(4.2) 



with 
dpo 



m^ + mg - 2M^ + 



(m^ 



rrin 



qi - 6momq2 



+ {n ■ no) 



m + itIq — 



(m — mo) 

w 



q2 + 2mmogi 



+ 2(no • p) [n ■ lo) q2 



+ 2(n • /o)t^(^' -ml- 2M^)g3 + 2(no ■ p)^{ml - - 2M^)g3 

My My 

In terms of the left and right-handed Ff^V couplings, the charges 51^2,3 are 



?1,3 



19± 
4 e2 



FfoV.2 



(12 



2 e2 



R 



(4.3) 



(4.4) 



Summing over the polarization of /o, the angular distribution dF/dcos^o where ^0 is the 
angle between the spin vector n and the flight direction of /o, can be written as 

dr 1 

-X- = 7^ (Ftot + cos^oF^g) (4.5) 

d cos 6*0 2 
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Ftot is obtained by integrating the previous expression over the angle 60 



tot 



d cos 6*0 — 

1 d cos 6*0 



(4.6) 



and corresponds to the partial decay width for on-shell gauge bosons. The cos^o term can 
be isolated by integrating dr/dcos6'o asymmetrically 



ang 



+1 dr 

dcosalo — 

d cos t/Q 



dr 

d cos 6*0 — 

-1 d cos 6^0 



The expressions of Ttot and Tang are 



Ttot 

r 

J- ang 



a m , 1 



2 M2 
a m? 
2M 



{ [(1 - f^iy + f^vi^ + /^o - 2/iy)] qi - 6/io/iyg2} 

■^(1 - /^O - 2/iy)g3 



where 



A = (1 - yUg - ;Uy) - AjjQjJy 



with /io = mo/m , /iy = My /m 



(4.7) 

(4.8) 
(4.9) 

(4.10) 



The partial decay width and the angular distribution in the case where the heavy fermion 
directly decays into its light partner and a gauge boson V can be obtained from the last 
expression by simply setting mo = in the previous expressions; one has 



Ttot 
r 

ang 



a Qy m? 



{Q 



R ) 



(4.11) 
(4.12) 



in agreement with Ref. ||2^, once the charges are specified. 



In the next subsection, we will discuss the case where the exchanged gauge bosons are off- 
shell, leading to three-body decays of the heavy fermions. 



4.2 Three Body Decays 



The amplitude for the decay, eq. (4.1), in the general case where the fermion /o is also massive 
and polarized can be obtained from the amplitude eq.(3.3) of the production of two heavy 
fermions by crossing symmetry: one simply has to change the labels of the four-momenta 
and to reverse the sign of the mass m. One obtains for the differential decay width 

dr = (2vr) 5{p-k-h- (2^)32/0 (2^)32/0 (2^)32/0 ^ (4.13) 
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where, after some simplifications, dr° reads 

2dr0 = {pd2){lo.k){Qi + Q3) + {p.h){lo.k){Qi-Q3)-rnmo{h.l2)Q2 

+ in.h)[moip.l2)Q'2 - m(/o./2)(Q'i - Qs)] - in.l2)[moip.li)Q'2 - m{k.h){Q\ + Q'^)] 

- {nQ.li)[m{lo.l2)Q'2 - mo{p.l2){Q'i + Q's)] + {nQ-h)[m{k.li)Q'2 - mo{p.li){Q[ - Q'^)] 
+ {n.l2){no.li)[m^Q2 + mmo{Qi + Q-^)] - {n.li){no.l2)[m'^Q2 + mmo{Qi - Q^)] 

+ {n.no)[{h.l2){p.lo) - (p./i)(/o./2) - {p.l2){lo.h)]Q2 

- 2[{n.h){no.l2) + in.l2)ino.l2mi.l2)Q2 (4.14) 

The generahzed charges Qi and Q[ are the same as those given for the production amphtudes, 
eq. (3.4), but they are built-up with the helicity amplitude 

assuming that the decay is mediated by only one gauge boson exchange; other channels can 
be easily included. 

In the case where the fermion /o is also massless, as it happens in most of the cases: the 
heavy fermion directly decays through mixing into its light partners and real or virtual 
gauge bosons which subsequently decay into two massless fermions, one can sum over its 
polarization and the expression dF'' simplifies considerably 

dr° = {p-l2){io-h){Qi + Q3) + {p-li){lo-h){Qi-Q3) 

-m{n ■ h)ilo ■ hW, - Q's) + min ■ /2)(/o ■ h)iQ[ + Q3) (4-16) 

However, since in this case several decay channels are possible, the helicity amplitudes are 
involved and one has in the general case where s = {p — Iq)"^ , t = {p — lif and u = {p — I2Y 
channels V5, Vt and Vjj exchanges are present^ 

2 /n/i/2Vs/-)/oFVs o 2 nfoh^T nf^^^T ^2 „2 nfofiVv nh^^u ^2 



„2 ^ _ T\/r2 „2 + _ T\,t2 



(4.17) 

where the widths of the gauge bosons have been omitted for simplicity. These expressions 
may need to be supplemented by statistical factors. 

However, in most physical situations there are at most two decay "channels" only, although 
in the same channel, several gauge bosons can be exchanged. These two decay channels 

■^An example of a situation where all the three channels occur is the decay of a heavy Majorana neutrino 
into an e^e~ pair and a neutrino or antineutrino. There is an "s-channel" Z-boson exhange N — > u^.Z — s- 
Vf.e'^e^ , a "t-channel" VF-boson exchange N — > e~W^ e+e^i/g and a "u-channel" exchange N 
e^W^ — )■ e^e~Ve] if the light neutrino is also a Majorana particle, the three amplitudes and the one due to 
N — > v^Z — > j/gC^e^ add coherently. 
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occur for instance, when two identitical particles are present in the final state [e.g. E — > 
e~Z — > e~e'^e~] hence requiring the symmetrization of the amplitudes or for decays involving 
both the isospin up and down light partners of the heavy fermions, leading to both neutral 
and charged gauge boson exchanges [e.g. E ^ e~ Z + z^eW^~ ~^ e~z/eZ7g]. In the rest of the 
discussion, we will therefore assume that only two channels are present and for simplicity, 
that there is only one gauge boson exchanged in each channel. 

Integrating over the momentum of one of the final fermions, e.g. /2, as well as on the 
azimuthal angle dependence, and using the usual scaled variables 



Xi 



2 (p • h) 



2 {p-k) 



the differential decay width writes 

dr 



dcos^idxidx2 



2 {p-lo) 



a 

Nc — m dTn 
Stt 



2 — Xi — X2 



(4.18) 



(4.19) 



with 01 the angle between the fermion /i and the spin vector n of the heavy fermion [the 
projections of n onto the fiight directions of /i and /2 are simply 2n ■ li/m — Xi cos^i and 
2n ■ 12/m = —X2 cos^i]; neglecting the widths of the exchanged gauge bosons and using the 
scaled masses vs — My^/w? and vt = My^/m^, dTg is given by 



xi{l - Xi)qZ^ + X2{1 - X2)q^ ^ Xo{l - Xo)q!J + ^2(1 - X2)ql 



TT 



(l-Xo- Vs)^ 



(1 



Xi 



2r]q^ 



{1-Xo- Vs)'^ 
ST ^2(1 - X2) 



+ 



{1-xi- Vt) 



[1 - xo - vs)(l - xi - Vt) 



cos 6*1 
(4.20) 



Here, is a statistical factor: 77 = —1 in the case where there are two identical Dirac fermions 
in the final state otherwise rj = +1; for quarks one has also to divide by the color factor, 
i.e. f] = ±1/3. The charges q^^ and q^'^ where M,N — S,T describe the channels with 
Vs, Vt exchange and the interference term; they read 



'MN 



1 r 

2 
1 

2 



M N 



M N 



111 QllI + Wrr Qrr 



\Qrr Qrr\ 





1 r 




~ 2 - 




1 r 




~ 2 



\1rl 1rl\ 



\Qrl Qrl\ 



IQlr 1lr\ 



\Qlr Qlr\ 



(4.21) 



where the helicity amplitudes qfj with i,j = L,R are given by 



_ 9Vs_^FfoVs^hf2Vs 



9Vt QFflVTQfof2VT 



(4.22) 
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Integrating over the energy of the particle /2 with the boundaries 1 — a^i < X2 < 1, one has 



dr 



dxidcos6'i 



{ x?(l - X]) 



^vr \vs{vs-X]) 



qSs _ q'ss 



5Xl/3) Irprp 

+ 7^ ^Mrrg cos^i 



+ 



+ 



(2x1 - 2vs - 1) log ^^^^ + ^(x,-l- 2vs) 

Vs Vs 

-2xi/3) 



qj' COS + qT_^ - qj^ cos 



2{1-Xi-vt)'^ 

- ^ cos 6i 
2 -qq^^ 



1 — Xi — Vt 



{1 - Xi + vs){vs - xi) log 



Vs 



- Xl[-Xi-1 - Vs 



(4.23) 



Finally, when the energy of /i is also integrated out, one has for the decay distribution 



dr 



d cos 9i 2 



- (Ftot + cos 6*11 



angj 



(4.24) 



where Ftot corresponds to the partial decay width for off-shell intermediate vector bosons 
and Tang to the angular distribution. They are given by 



Ttot 

r 

ang 



8n 



AT 

— l\r ^ 

"8n 



{q^J + gf ) Rt{vs) + (g^^^ + qr) RtM + Vql^Riivs, Vt)] (4.25) 
(gl^^ + q'f) Rt{vs) + (gl^^ + q'r) RtM + qTRAM] (4.26) 



where Rt{v) describes the decay width in a given channel, 

V — 1 1 

R^fv) = v(v - 1) log + — (6t;^ - 3t; - 1) 

V 6f 



(4.27) 



Ri{vs, Vt) describes the interference between the two channels when simultaneoulsy present. 
In terms of the Spence funtion, Li2(x) = — Jq dy log(l — xy), it is given by 



Ri{v,v) 



[1 - VT)il - 2vs - ^VT)\og 



Vt 



Vt — I 



+ {l-vs){l-2vT-3vs) 



Vs 

X log - + 2{vs + VT)il-Vs-VT) 

vs-l 



Vs + Vt - I 



Vt + vs - I 



■ Vt , vs + vt-1 
log log 

Vt — I Vs 

+ 3 — 5^5 — 5vt 



(4.28) 



and Ra describes, for only Vt exchange, the deviation of the angular distribution from the 
familiar dr/dcos6'i ~ (1 — cos 6^1) form. 



RAiv) = - 1)(3 - 5v) log ^ - -^(6t;2 - 21v - 8) 
2 V 12v 



(4.29) 



Once the charges are specified, the expressions of Ry and Rj agree with those obtained in 



Ref. [|T9| and the expression of Ra [and Ry] with the result of Ref. pO 
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5. Summary 



In this paper, we have analyzed in detail the decay modes and the production mechanisms 
in e"'"e~ collisions, of new heavy fermions predicted by extensions of the Standard Model and 
in particular by Grand Unified Theories. 

We have given analytical expressions for the production of two heavy fermions of any flavor 
in a general case: we have allowed for arbitrary couplings of the new fermions, the presence 
of several production channels and we have taken into account the polarization of the initial 
e^e~ state and the final polarization of the heavy fermions. We have treated the case 
where the two fermions have different masses, to discuss the pair production and the single 
production in association with ordinary fermions on the same footing, and to account for 
the possibility of producing two different heavy leptons. Complete and compact formulae 
were given for angular distributions, total cross sections and polarization vectors. 

We have then discussed in detail the decay modes of the new fermions including the cas- 
cade decays and the three-body decays with off-shell gauge bosons. Complete analytical 
expressions were given for total widths as well as for angular and energy distributions. We 
have also taken into account the polarization of the decaying fermion; combined with the 
spin-dependent cross sections, the spin-dependent decay distributions allow to obtain the 
full correlations between all the particles involved in the process. These correlations are very 
useful to discriminate between different types of particles and will be very important when 
discussing the signals for heavy fermion production since they help to optimize the experi- 
mental cuts which permit to suppress the various backgrounds without affecting drastically 
the signal cross sections. 

More phenomenological aspects of heavy fermion production at future high-energy e"'"e~ 
linear coUiders, including a detailed analysis of the various signals and backgrounds, will be 
presented in a subsequent paper p3| . 
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Appendix: Correlated Production and Decay 



In this Appendix, we briefly summarize the formahsm introduced by Tsai to describe 
the production of heavy fermions in e~^e~ coUisions 

e+e" — > F F (A.l) 

which subsequently decay through real or virtual gauge boson exchange, into three body 
massless final states 

F — ^ /o /l /2 

F ^ Tohf2 (A.2) 

including the spin correlations. This formalism allows an easy reconstruction of the correla- 
tion between the initial e"*", e~ states and the final decay products of the heavy fermions. 

The polarization vector V{V) of the fermion F{F) is defined with the help its spin four- 
vector n{n). The latter satisfies the relations n ■ n = n ■ n = —1 and n ■ p = n -p = 0, and is 
introduced by replacing the usual projection operators by 

u{p)u{p) {^ + m) ^ for particles 
1 -|- 75"^ 

v{p)v{p) —>■ — m) for antiparticles (A. 3) 

In the production process, the polarization vector V is defined as 

d(TP°^ = ^ da'^^P"' + (A.4) 

Choosing the (x, z) plane as the scattering plane with the electron along the +z direction, 
the covariant spin vector n can be decomposed, in general, along three directions: the F 
direction p {n\\), the transverse direction with respect to p but within the scattering plane 
(n1) and the transverse direction with respect to p but along a normal to the scattering plane 
{um)- The projection of the spin vector along these three directions defines the corresponding 
degrees of polarization 

= i^infl + P^< + PNn% (A.5) 

In the case where there is no CP violation, as we will assume for the process (Al), and since 
the very small imaginary parts from width or loop effects can be safely neglected, there is no 
polarization transverse to the scattering plane, Pn = 0. We can therefore set the azimuthal 
angle to zero and take = {E, \p\ sin 9, 0, |p| cos 6') where 9 is the scattering angle; in this 
case the components of the polarization vector are simply 

V\\ = E/m {\p\/E,sm9,0,cos9) , = (0, cos^, 0, - sin ^) (A.6) 
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Taking into account the polarization of both final fermions, the cross section for the pair 
production, eq. (Al), can be written as 

dfjP"^ = ^ d(7""P (l + V-n + V-n + C^'^n^n,) (A.7) 
with C^^ the spin correlation. 

In the decay processes, eq. (A2), the polarization vectors of the heavy fermions V', are 
defined similarly to eq. (A2) but without the factor 1/2; for instance 

(ippol ^ ^punpol (1 + ^ . 7?') (A.8) 

In the narrow width approximation, the full cross section for the production and subsequent 
decay of the heavy fermions, 

e+e-^/o^/i (A.9) 

can be written as 

da = id(j""P°'dr"'^P°' dr"P°' (l + rii,, V^V" + rJ^JP^W + V^aVr^^d^'V'"!^^) (A.IO) 

with 

ViMv = -Qixu +PnPu/m^ , r)^^ = -Q^iv +P„P^/rff (A.ll) 
Summing over the polarization of one fermion, e.g. F, the full cross section simphfies to 

d(7 = ^ d(7"°P°' dr""P°^ (l + r]f,,V^V'") (A.12) 

— * — * 

Note that iJni^V^V"^ — V*.V'*, where the * refers to the components in the rest frame of the 
heavy fermion. 

The above formalism, based on the factorization of the production and decay sequences, 
permits an easy reconstruction of the full correlations between the initial and the final 
particles from the decay of the heavy fermions, which would be very difficult to obtain 
directly from the full process eq. (A9) with six particles in the final state. Furthermore, it 
is well adapted for the setting of a Monte Carlo generator. 
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